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ABSTRACT 

Aims. We investigate the X-ray spectral and timing properties of the high mass X-ray binary EXO 2030+375 observed during its 
lune-September 2006 giant (type II) outburst. 

Methods. The data analyzed in this work are from partly simultaneous observations with INTEGRAL and Swift. The pulse period P 
and its temporal derivative P are measured. X-ray pulse profiles in different energy ranges and time intervals are constructed. Pulse 
averaged X-ray spectra for different time intervals are studied. 

Results. We report a strong spin-up of the source during the outburst, comparable to that observed in 1985 during the previous giant 
outburst when the source was discovered. The value of P is found to be linearly related to the X-ray luminosity of the source during 
the outburst. For the first time the hard X-ray (> 25 keV) characteristics of the source during a type II outburst are studied. The X-ray 
pulse profiles apparently change with luminosity. The X-ray spectral continuum in the 3-120 keV energy range is modeled with an 
absorbed power law with an exponential cutoff around E ~ 26 keV. An iron emission line at ~6-7 keV is observed. The spectrum 
reveals some features between 10 and 20 keV which can be modeled either by a broad emission line at ~13— 15 keV (a "bump") or by 
two absorption lines at ~10 and ~20 keV. 
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> 1. Introduction 

• i-H , 

| Discovered by EXO SAT in 1985 during a giant outburst 
■ dParmar et al.l 11989b . the transient accreting pulsar EXO 
2030+375 belongs to a high mass X-ray binary sys tem with 
a B0 Ve star as optical companion dCoe et alj 1 1988b . The or- 
bital period and ecc entricity are ~46 d and ~ 0.42, respectively 
dWilson et al.ll2.Q05b . According to the generally accepted model 
for Be/X-ray binaries, the X-ray emission in such systems dur- 
ing outbursts is produced when the compact object accretes from 
a quasi-Keplerian disk around the equator of the rapidly rotating 
Be star. Such a mechanism explains normal (type I) outbursts 
with X-ray luminosity Lx < 10 37 ergs~', associated with the 
periastron passage of the compact object. Sometimes, however, 
Be/X-ray binaries show giant (or type II) outbursts with X-ray 
luminosity Lx > 10 37 ergs~'. Such outbursts are thought to be 
due to a dramatic expansion of the disk surrounding the Be star, 
leadin g to the form ation of an accretion disk around the compact 
object dCodl2000l) . 



Table 1. Summary of observations. 



Obs. 


Instrument 


Obs. time 
(ks) 


Mean 
MJD 


Mean ASM 
flux (mCrab) 


OBS1 


ISGRI+JEM-X 


62 


53942.9 


500 


OBS2 


ISGRI+JEM-X 


140 


53967.6 


500 


OBS2 


XRT + BAT 


6.3 


53967.3 


500 


OBS3 


XRT + BAT 


6.4 


54002.5 


160 



* Based on observations with INTEGRAL, an ESA project with in- 
struments and science data centre funded by ESA member states (espe- 
cially the PI countries: Denmark, France, Germany, Italy, Switzerland, 
Spain), Czech Republic and Poland, and with the participation of Russia 
and the USA. 



During the 1985 giant outburst the X-ray luminosity of the 
source reached a value of £1-70 fev ~ 2 • 10 38 erg s (for a dis- 
tance of 7.1 kpc, Wils on et all 12002). The spin period of the 
pulsar (~ 42 s) changed dramatic ally, with a spin-up time scale 
-P/P * 30 yr dParmaretalJll989l) . 

From the nearly continuous monitoring of the source with 
the Burst and Transient Source Experiment (BATSE) on the 
Compton Gamma Ray Observatory (CGRO) and the All Sky 
Monitor (ASM) on the Rossi X-ray Timing Explorer (RXTE) 
it is known that the pulse period remained roughly constant for 
about a year after the beginning of the monitoring in 1991, fol- 
lowed by 2 years of relatively slow s pin-up and 6 years of slow 
spin-down dWilson et al .1120021 [2005b . After 2002 a transition to 
a global spin-up along wi th an overall brighten ing of the normal 
outbursts was reported bv lWilson et al.1 d2005l) . 
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Fig.l. The ASM light curve of EXO 2030+375. The times 
of INTEGRAL (INT.) and Swift observations analyzed in this 
work are marked by solid lines. The dashed line shows another 
simultaneous observation by INTEGRAL and Swift for which 
the data are not yet available. 



The giant outburst of EXO 2030+3 7 5 presented here 
dCorbet & Levind 120061: iKrimm et al] 120061: iMcCollough et al.1 
2006) is the second of its kind since the discovery of the source. 
The X-ray luminosity at the m aximum of the outbur st reached 
L1-20 kev ~ 1.2 • 10 38 ergs' 1 dWilson & Finger! |2006|) . In this 
Letter, based on INTEGRAL and Swift data, we investigate in 
§ 2.1 the pulse period and the spin-up rate of the source and 
in § 2.2 we construct X-ray pulse profiles for different energy 
ranges and time intervals. To characterize the broad-band spec- 
tral behaviour we perform in § 2.3 the spectral analysis of pulse- 
phase averaged spectra of the source. A detailed study of the 
timing and spectral behaviour as a function of the luminosity is 
beyond the scope of this Letter and will be presented in a forth- 
coming paper. 



2. Observations and data analysis 



EXO 2030+375 was observed by INTEGRAL dWinkler et alJ 
12003b on 27 July and 19-21 August 2006 and by Swift 



(Geh rels etalll2004|) on 19-20 August and 23-25 September 
2006. Another simultaneous observation by INTEGRAL and 
Swift was performed on 6-8 October 2006, but the data were not 
available at the time of writing this Letter and will be presented 
later. For our analysis we used the data obtained with the instru- 
ments IBIS/IS GRI (20-300 keVjUbertini et al.l l2003l) and JEM- 
X (3-30 keV. ILund et alj|2003l) onboard INTEGRAL as well as 
XRT (0.2-10key.lBurrows et al.l l2005l) and BAT (15-150 keV, 
Bart helmv et al.ll2005l) of Swift. For the analysis of ISGRI and 
JEM-X data the Off- line Science Analysis (O SA) software (ver- 
sion 5.1) was used (Co urvoisier et al.ll2003l) . Analysis of XRT 
data was done with the standard HEADAS package (version 
6.0.4). 

Table Q] contains the summary of the observations analyzed 
in this work. A part of the ASM light curv <0 of the source in- 
cluding three normal outbursts and the 2006 giant outburst is 
presented in Figure Q] The observations analyzed in this paper 
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Fig. 2. Spin frequency of EXO 2030+375. Small diamonds show 
the spin frequency de velopment before the 200 6 giant outburst 
(they are taken from IWilson et al.ll2002l 12005). The large tri- 
angles correspond to our pulse frequency measurements during 
the outburst. The uncertainties of measurements and duration of 
each observation are smaller than the size of symbols in the plot. 

Table 2. Pulse period development. The uncertainties in paren- 
theses (68%) refer to the last digit(s). 



Obs. 


MJD 


P 
[s] 


-dP/dt 
[10- 8 s/s] 


Flux (3-10 keV) 
[10~ 9 ergcirT 2 s -1 ] 


OBS1 


53942.7714 


41.57958(2) 


2.93(12) 


9.61(6) 


OBS2 


53967.4260 


41.51706(5) 


3.29(19) 


10.41(4) 


OBS3 


54002.3877 


41.45954(4) 


0.67(10) 


3.25(3) 



1 We used the results provided by the ASM/RXTE teams at MIT and 
at the RXTE SOF and GOF at NASA's GSFC. 



are marked by solid lines and consist of three sets of pointings. 
The first set was done by INTEGRAL. The second set contains 
simultaneous observations by Swift and INTEGRAL. The last 
set of observations was performed by Swift. These three sets will 
be referred to as OBS1, OBS2, and OBS3 throughout the paper. 
OBS1 and OBS2 observations were made when the source was 
at approximately the same luminosity level, before and after the 
maximum of the outburst. OBS3 set was made during the decay 
of the outburst when the luminosity of the source dropped by a 
factor of ~3. 



2. 1 . Pulse period behaviour 

For the timing analysis all times were translated to the solar sys- 
tem barycenter and corrected for orbital motion in the binary. 
The pulse periods and associated derivatives were determined 
for each of the three observation periods individually (Tabled. 
These values were found by employing initial epo ch-folding and 
a subs equent phase connection analysis similar to lFerrigno et all 
(2006) using well defined pulse profiles from a sufficiently large 
number of pulses. Variation in pulse shape inside each observa- 
tion is marginal and does not affect our method. 

Fig. 2 displays our points toget her with the h i storical puls e 
period development as presented by Wilso rTet alJ d2002, 2005). 
The three datasets are consistent with a global solution (valid 
between MJD 53943 and MJD 54003) with the following pa- 
rameters: P = 41.47268(6) s, P = -1.627(5) ■ 10~ 8 s/s and 
P = 4.058(25) ■ 10~ 15 s 1 for the reference MJD 53992.25 
(TDB). The uncertainties in parentheses (68%) refer to the last 
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Fig. 3. X-ray pulse profiles obtained with ISGRI (top panel) and 
JEM-X (bottom panel) during OBS2 (MJD- 53966), near the 
maximum of the outburst. Phase 0.0 was arbitrarily chosen to 
be defined by the dip-like feature before the main peak. Typical 
error bars are indicated. 

digit(s). The characteristic spin-up time scale, -P/P, is ~40 
years, somewhat longer than during the first giant outburst. We 
also note that the P appears to be linearly correlated with the 
X-ray luminosity (see last two columns in Table|2]i. The formal 
correlation coefficient is 1.0 and the slope from a linear fit is 
36 + 2 (s s~')/(erg cm _2 s _1 ), although it is based on three data 
points. 

2.2. Pulse profiles 

Using the values of P and P determined here (Table |2| we con- 
structed pulse profiles for each observation. Pulse profiles from 
OBS1 and OBS2 are found to be very similar. This can be ex- 
pected considering that the luminosity of the source during these 
observations was approximately the same. The results are shown 
in Figs. [3] and |U We do not show pulse profiles from OBS1 since 
there are similar to OBS2. 

Fig. |4] shows the 0.2-10 keV pulse profiles obtained with 
XRT. One can see that the pulse profile obtained in OBS3, dur- 
ing the decay of the outburst, appears different from that ob- 
tained in OBS2 when the luminosity was higher. 

The analysis of the pulse profiles, including the comparison 
with those obtained during the 1985 outburst and during normal 
outbursts, is on-going. 

2.3. Spectral analysis 

For OBS 1 pulse averaged X-ray spectra of ISGRI and JEM-X 
were analyzed. For OBS2 spectra from ISGRI, JEM-X, XRT 
and BAT were fitted simultaneously. In order to account for 
small scale uncertainties in the response matrices of the respec- 
tive instruments systematic uncertainties were added at a level of 
1% for ISGRI and 3% for JEM-X and XRT. For BAT we used 
energy-dependent systematic errors provided by the HEADAS 
package. We note that there are additional global systematic un- 
certainties with the instruments on INTEGRAL (e.g. the canon- 
ical spectrum of the Crab is not reproduced). The stated uncer- 
tainties (e. g. in Tables 3 and 4) are therefore to be taken as lower 
limits. 

To fit the broad band continuum of the source we used an ab- 
sorbed power law with an exponential cutoff (XSPEC highecut 



model): 

dN^fE' 1 ', if£<£cutoff 
dE K \E r ■ exp if £ > £ cutoff , 

v ^folil 

where dN/dE is the differential count rate, E is the energy of 
the photons; T, E cuto ff, and Ef \d are model parameters. An iron 
emission line at ~6.5 keV was observed during OBS2 (which has 
the largest total number of counts) and modeled by a Gaussian 
emission line. JEM-X data revealed an additional feature in 10- 
20 keV range, which could be explained by a Gaussian emission 
model at ~ 1 3—15 keV (a "bump"). Alternativ ely, following the 
sugges tion of a cyclotron line at ~10 keV by IWilson & Finger! 
(2006), we equally well fitted the spectra adding two Gaussian 
absorption lines at ~10 and ~20 keV to the broad band contin- 
uum. To account for large systematic uncertainties in the abso- 
lute flux measured by the instruments we introduced in our mod- 
els a free multiplicative factor for each instrument: FisGRb ^xrt> 
and Fbat (for JEM-X the factor was fixed to 1.0). The best fit 
spectral parameters for OBS 1 and OBS2 are listed in Table[3](for 
the model with a "bump") and Table [4] (for the model with two 
absorption lines). The broad band spectrum of OBS2 fitted with 
highecut is shown in Fig.|5^. Figures^, HJ; and [3Jl show the 
residuals after fitting the spectrum by the highecut model with- 
out additional features between 10 and 20 keV, with a "bump" 
around 15keV, and, alternatively, with two Gaussian absorption 
lines. 



3. Summary and Discussion 

In June-September 2006 EXO 2030+375 has undergone the sec- 
ond giant (type II) outburst since its discovery in 1985. During 
the outburst the source has been observed with INTEGRAL and 
Swift. The main results of these observations can be summarized 
as follows: 

1. The spin frequency of the pulsar has dramatically in- 
creased, indicating the presence of an accretion disk around the 
neutron star. The maximum value of P is found to be compara- 
ble to that observed during the first outburst during which the 
source was discovered. The values of P determined for the three 

Table 3. Best fit spectral parameters of EXO 2030+375 using 
the model with a "bump" around 15 keV (see text). l<x(68%)- 
uncertainties for one parameter of interest are shown. 



53943(OBSl) 53967(OBS2) 



Time of observation 


ISGRI, 


ISGRI, IEM-X 


(MJD) / Instruments 


JEM-X 


XRT,BAT 


r 


1.93 ±0.01 


1 QQ+0.01 

1-7 -0.02 


£cutoff [keV] 


25 9 +0 - 2 

-0.3 


26 4 +0 2 

iu - -0.5 


£foid [keV] 


26.1 ±0.2 


26.9 ± 0.2 


£ Fc [keV] 




6.6 ± 0.1 


o-Fe [keV] 




1.2 ±0.1 


Ebump [keV] 


15.3 ±0.2 


13 4 +04 


0-bump [keV] 


2.7 ± 0.2 


4 l +0 3 


N K [10 22 cnT 2 ] 


3.1 ±0.2 


3.4 ±0.1 


F JEM-X 


1.0 (fixed) 


1.0 (fixed) 


F ISGRI 


1.20 ±0.01 


1.20 ±0.01 


^XRT 




0.926 ± 0.004 


^BAT 




0.91 ±0.01 


xlJA.o.f. 


1.1/259 


1.3/937 
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Fig. 4. X-ray pulse profile obtained with XRT during OBS2 (top 
panel) and during OBS3 (bottom panel). Typical error bars are 
indicated. 



observation periods are linearly rel ated to the X-ray l uminosity, 
in line with previous observations (Parmaret al. 1989) and the 
prediction of accretion torque theory dGhosh & Lambl[l979l) . 

Since the luminosity of the source and its spin-up rate could 
be measured, we can estimate the value of the magnetic field 
strength on the neutr on star using the clas sical model of an ac- 
creting X-ray pulsar dPringle & Reesl[l972l) . The equation deter- 
mining the change of the angular momentum of the neutron star 
has the form Ico — M yjGMR^, where / is the moment of in- 
ertia of the neutron star, a> is the angular frequency, Ra is the 
Alfven radius, M is the mass of the neutron star, and M is the 
accretion rate. Ra depends on the magnetic field strength and 
the accretion rate as Ra - [p 2 /(2M V2GM)] 2/7 , where ,u is the 
magnetic dipole moment of the neutron star. Using the observed 
values of the luminosity L (which can be translated into the ac- 
cretion rate using the relation L ~ 0.1 Mc 2 ), P and standard neu- 
tron star parameters, we obtain an estimate of the magnetic field 



Table 4. Best fit spectral parameters of EXO 2030+375 using the 
model with two absorption lines at ~10 and ~20 keV (see text). 
lcr(68%)-uncertainties for one parameter of interest are shown. 





53943(0BS1) 


53967(OBS2) 


Time of observation 


ISGRI, 


ISGRI, JEM-X 


(MJD) / Instruments 


JEM-X 


XRT,BAT 


r 


1.79 + 0.01 


1.78 + 0.01 


Ecutoff [keV] 


19.8 ± 0.4 


19.5 + 0.3 


£foid [keV] 


24.2 ± 0.2 


24 7 +al 

^•'-0.2 


£ Fe [keV] 




6.51+0.03 


O-Fe [keV] 




n ^9+0.04 
"- 3i _0.03 


Eist line [keV] 


io.o!° : 2 


10 6 +al 


0-isi lmc [keV] 


1-6^ 


0.7 + 0.2 


E 2 na line [keV] 


20.7 ± 0.3 


20.6 ± 0.3 


0" 2n d line [keV] 


1.3 + 0.3 


2 +03 


N H [10 22 cnr 2 ] 


2.2 + 0.2 


3.1 + 0.1 




1.1/256 


1.3/934 


^JEM-X 


1.0 (fixed) 


1.0 (fixed) 


flSGRI 


1.25 ± 0.01 


1.20 + 0.01 


^XRT 




0.926 ± 0.004 






0.90 ± 0.01 






10 20 
channel energy (keV) 

Fig. 5. The broad band spectrum of EXO 2030+375 from si- 
multaneous fits of INTEGRAL and Swift data of OBS2 with 
highecut (a) and residual plots after fitting it without additional 
features (b), adding a "bump" around 15 keV (c), or alternatively 
including two absorption lines at ~10 and ~20 keV (d). 



strength B ~ 1 - 4 • 10 12 Gauss. This would be consistent with 



a fundamental cyclotron line in the range E, 



eye 



10 - 50 keV 

(E cyc = 11.6 ■ (B/10 12 G) x [1 + zr 1 keV, where z is the gravita- 
tional redshift which was assumed to be 0.2 in our calculations). 

2. For the first time the broad band X-ray spectrum of the 
source (including energies above ~25 keV) during a type II out- 
burst is analyzed. The continuum is modeled by a power law 
with an exponential cutoff. An iron emission line is observed at 
~6-7 keV. The spectrum shows additional features between 10 
and 20 keV which can be equally well modeled with a broad 
emission line at ~ 13 - 15keV or two absorption lines at ~10 
and ~20 keV. On the basis of phase averaged analysis we can not 
distinguish between the two options and we can not c onfirm the 
presen ce of a cyclotron line at ~36 keV reported by Reig & Cod 
( 1999). A preliminary pulse phase resolved analysis shows that 
the spectrum is highly variable with pulse phase. This means 
that any interpretation of the phase averaged spectrum should be 
made with caution. Moreover, residual features around 25 keV 
may hint at not yet well understood calibration uncertainties. 

3. The shape of the profiles is found to vary during the out- 
burst, probably related to the luminosity. 

Pulse-phase-resolved analysis along with the analysis of the 
spectral and timing behaviour of the source at different luminos- 
ity states, is on-going. 
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